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We provide conservative bounds on the dark matter cross-section and lifetime from final state 
radiation produced by annihilation or decay into charged leptons, either directly or via an interme- 
diate particle (j>. Our analysis utilizes the experimental gamma-ray flux upper limits from four Milky 
Way dwarf satellites: HESS observations of Sagittarius and VERITAS observations of Draco, Ursa 
Minor, and Willman 1. Using 90% confidence level lower limits on the integrals over the dark matter 
distributions, we find that these constraints are largely unable to rule out dark matter annihilations 
or decays as an explanation of the PAMELA and ATIC/PPB-BETS excesses. However, if there 
is an additional Sommerfeld enhancement in dwarfs, which have a velocity dispersion ~ 10 to 20 
times lower than that of the local Galactic halo, then the cross-sections for dark matter annihilating 
through <j}'s required to explain the excesses are very close to the cross-section upper bounds from 
Willman 1. Dark matter annihilation directly into r's is also marginally ruled out by Willman 1 
as an explanation of the excesses, and the required cross-section is only a factor of a few below 
the upper bound from Draco. Finally, we make predictions for the gamma-ray flux expected from 
the dwarf galaxy Segue 1 for the Fermi Gamma-ray Space Telescope. We find that for a sizeable 
fraction of the parameter space in which dark matter annihilation into charged leptons explains the 
PAMELA excess, Fermi has good prospects for detecting a gamma-ray signal from Segue 1 after 
one year of observation. 



I. INTRODUCTION 

While the evidence for non-baryonic dark matter is 
compelling on astrophysical scales, its microscopic nature 
remains unknown. One appealing possibility is that dark 
matter consists of particles interacting via the weak force 
(WIMPs). Such particles have an annihilation cross- 
section at freeze-out of ~ 3 x 10^^^ cm'^s^^, and if these 
particles are in the mass range of about 10 GeV to 10 
TeV, then this would yield the correct relic abundance of 
dark matter observed today. In addition, the hierarchy 
problem of particle physics suggests new physics beyond 
the Standard Model should enter on scales of 100 GcV 
to 1 TeV, providing additional motivation for the WIMP 
hypothesis. If dark matter consists of WIMPs, then it 
could annihilate or possibly decay into lighter particles 
which could be observed as gamma-rays, cosmic rays, 
or neutrinos. There are a host of current and planned 
instruments that have the potential to detect this sig- 
nal via an observed excess of these annihilation or decay 
products. Such indirect search techniques complement 
direct dark matter detection experiments and attempts 
at accelerator-driven dark matter production. 

Recent results from the PAMELA satellite indicate an 
excess of positrons between 10 GeV and 80 GeV but 
no excess in anti-protons 0, above the background ex- 
pected from cosmic ray nuclei interactions with the inter- 
stellar medium. This confirms earlier hints of a positron 
excess seen by other experiments (e.g., 0,13, SB)- The 
balloon-borne experiment ATIC recently reported an ex- 
cess of cosmic ray electrons plus positrons in the en- 
ergy range of 300-800 GeV [31, confirming earlier ex- 



cesses seen by PPB-BETS 3- Since cosmic ray elec- 
trons and positrons loose energy rapidly through in- 
verse Compton and synchrotron processes, the source of 
these excesses should be less than 1 kpc away Q. It 
is possible for these excesses to be explained by con- 
ventional astrophysical objects such as nearby pulsars 
[13, [O, [H, [H, [H- It is also possible that they may be 
the first indirect evidence of dark matter annihilations 
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In light of the latter possibility, it has been shown that 
both PAMELA and ATIC/PPB-BETS, as well as a pos- 
sible excess in microwave photons at the Galactic center 
dubbed the WMAP haze [H, [H, [13, [11, may be ex- 
plained by a ^ 1 TeV dark matter particle which anni- 
hilates into charged leptons [l^, [lB|- However, for this 
scenario to work, the annihilation cross-section in the 
Galactic halo must be boosted by factors of O(IOO) be- 
yond the freeze-out cross-section. One suggestion is that 
this could be achieved via a velocity dependent cross- 
section throu gh t he mechanism of Sommerfeld enhance- 
ment [l^, [TtI. ISGI. [stI . [sst . This idea supposes that the 
0(100) lower velocity dispersion in the Galactic halo as 
compared to the velocity at freeze-out could account for 
the necessary boost to explain these recent observations. 

The main signal of interest in this work is the gamma- 
ray flux from final state radiation (FSR) by charged lep- 
tons produced in dark matter annihilations or decays. 
Here FSR is defined as photons directly radiated from the 
external legs (charged leptons) of the interaction. This 
signal has the advantage that the gamma-ray spectrum 
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has a sharp cutoff at energies equal to the dark matter 
mass and tliat the normahzation and shape of the spec- 
trum can be predicted without knowledge of the astro- 
physical environment [soj . In addition, for dark matter 
masses of a few hundred GeV, the FSR signal is relevant 
at high photon energies (> 100 GeV) making it accessible 
to Atmospheric Cherenkov Telescopes (ACTs). 

We focus on dwarf galaxies because the largest FSR 
signal is expected from regions with the highest dark mat- 
ter density. The Galactic center is a promising source, 
however it is difficult to separate dark matter produced 
gamma-rays from other astrophysical backgrounds to- 
ward the Galactic center, and there is wide uncertainty 
in the shape of the central dark matter density profile. 
However, dwarf galaxy satellites of the Milky Way are 
relatively nearby, are highly dark matter dominated with 
mass to light ratios of 0(10 — 100), and are largely free 
of astrophysical backgrounds as they have little warm or 
hot gas, minimal dust, and no magnetic fields [60| . In 
addition, their dark matter distributions can be inferred 
directly from stellar kinematics, and their velocity dis- 
persions are of 0(10) less than in the Galactic halo (see 
e.g. [6l[), the latter of which could result in an order of 
magnitude larger boost of the annihilation cross-section 
if Sommcrfeld enhancement is relevant. 

In Section |TT1 we discuss the calculation of the gamma- 
ray signal from dwarf galaxies and provide the formulae 
we use to calculate the FSR signal. Section IIIII reviews 
the most recent gamma-ray constraints from the ACTs 
HESS and VERITAS. In Sections |TV] and El we give con- 
servative bounds on the dark matter annihilation cross- 
sections and decay lifetimes, respectively, for the channels 
in which dark matter goes to e'^e~ or n'^fJ-~ either di- 
rectly, or through an intermediate particle (resonance), 
which we call cj). We also consider the case in which 
dark matter goes directly to t~^t~. These bounds are 
implied by gamma-ray observations of the dwarf galax- 
ies Sagittarius, Willman 1, Draco, and Ursa Minor. We 
use updated stellar kinematic data to calculate the dark 
matter density profiles of these dwarf galaxies, marginal- 
izing over a wide range of profile shapes. We also dis- 
cuss the implications of our constraints for the dark mat- 
ter explanation of the PAMELA and ATIC/PPB-BETS 
excesses. In Section IVIl we discuss the prospects for 
the Fermi gamma-ray space telescope to detect a dark- 
matter-induced gamma-ray signal from the dwarf galaxy 
Segue 1. Among the known dwarf galaxies. Segue 1 could 
potentially provide one of the largest signals from dark 
matter annihilations (6^ . We summarize our results in 
Section Eni 



II. GAMMA-RAYS FROM ANNIHILATING 
AND DECAYING DARK MATTER 

The gamma-ray flux from annihilating dark matter in 
a dark matter halo is given by 
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where is the dark matter particle mass, {av) is the an- 
nihilation cross-section, Eth is the threshold energy of a 
given gamma-ray instrument, and E^ax is the maximum 
energy of the photons. The integral over depends 
only on the particle physics details of the dark matter 
annihilation. £ann is given by 
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and depends solely on the properties of the dark mat- 
ter halo and the solid angle over which it is observed. 
Here p{r) is the halo dark matter density profile, r = 



/D2 



— 2sDcos9 for a halo at a distance D, and the 



integral is performed along the line of sight over a solid 
angle Afi = 2tt{1 — cos 6*). 

We will also be considering constraints on decaying 
dark matter from FSR. In this case, the 7-ray flux is 
analogous to the annihilating case with the following sub- 
stitutions: {av)/{2M^) T/M^ and in £, ^ p. We 
define the equivalent integral to eq. ([2]) over the density 
for the case of decays as >Cdcc- Since decaying dark matter 
is sensitive to p and not p^, the dark matter constraints 
from gamma-ray observations will be weaker for the case 
of decays. 



A. Dwarf Galaxy Dark Matter Density Integrals 

We determine the line-of-sight integral in eq. ([2]) from 
the kinematic data in the dwarf galaxies, in a manner 
similar to the analysis of [1^ . For each galaxy we obtain 
the distribution of the factors >Cann and £dcc assuming a 
five-parameter density profile of the form 
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We marginalize over a,b,c, pQ,rQ, with uniform priors, 
over conservative and plausible values for these param- 
eters. The specific ranges for the shape parameters are 
a = [0-1.5], b = [0.5-1.5], c = [2-5]. Though numerical 
simulations produce cuspy profiles for cold dark matter 
(CDM) models, which profiles have a ~ 1 [g^], we allow 
for a broader range to account for the fact that the dark 
matter we consider here may be produced non-thcrmally, 
may be self-interacting, or may be warmer than standard 
CDM. In this way, our analysis differs from [6^, in which 
CDM priors and cuspy profiles were assumed. For Sagit- 
tarius, Ursa Minor, and Draco, the scale radius is fixed 
to be within the range tq = [0.1 — 10] kpc, while for Will- 
man 1 the scale radius is fixed to be within the range 
To — [0.01 — 1] kpc. For each galaxy, po is set by the 
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condition that the mass within 300 pc is fixed to the re- 
spective observed values (65| . It is important to note that 
even though we assume a wide range for the parameters 
that determine the density profile, the velocity dispersion 
data fix the distribution of the line-of-sight integral to lie 
within a relatively narrow range. This is true despite 
the fact that the parameters (a, 6, c, po? ^o) SlTC not well 
constrained in and of themselves. This fact precludes us 
from showing the best-fit values for {a,b, c, pq^tq). Re- 
garding the difference between cored and cupsed models, 
we find that allowing for cored models only reduces the 
90% confidence level lower limits, which we use below, 
by ~ 10% for each galaxy relative to a ~ 1 models. This 
traces back to the fact that the stellar kinematics fixes 
the normalization of the density and density-squared for 
the dwarf dark matter halos, so that varying the inner 
slope has little impact on the final results p3^ . 

The 1-dimensional stellar velocity dispersions, aver- 
aged over each of the respective galaxies, are listed in 
Table ID We find that the velocity dispersion of the stars 
is similar to that of the dark matter for these dwarf galax- 
ies. Assuming the dark matter velocity dispersions are 
isotropic, we estimate the 3-dimensional dark matter ve- 
locity dispersions to be -\/3 times these values. 

Two galaxies listed in Table [I] are worth commenting 
on. First, the satellite Willman 1 was discovered in 2004, 
and while there is evidence that this object is dominated 
by dark matter, better kinematic measurements are re- 
quired to provide constraints on its dark matter distribu- 
tion [gI]. Second; it is clear that Sagittarius is currently 
being tidally-disrupted, and for an equilibrium model to 
be self-consistent the stars used in the analysis must be 
bound to the central core of the galaxy. To conserva- 
tively account for this issue, in our analysis we only use 
the stars from Ref. [131 within the bound region of ~ 0.4 



kpc [6^. The remaining two dwarfs. Ursa Minor and 
Draco, are conclusively known to be dark matter domi- 
nated E^M. 



B. Final State Radiation Calculations 

For FSR, for the dark matter annihilation channel 
XX ^ ff^ where / is a fermion, the photon energy spec- 
trum is given by 
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Here a ~ 1/137, y = E^/m^, 
= in eq. (HI). This formula 
holds in the coUinear limit, where the photon is emitted 
coUinearly with one of the leptons and m/ <C m^. For 
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the decay channel x 



2E^/m 



ff, we have instead y 
s = m^, and E^ax = "i^/S. 

To_ calculate the FSR for the channel XX ~^ 00 ^ 
//// [H, [m, [ill , where (/) is a light particle (resonance) 



we start with eq. (j4|) . where now s 



and y is replaced 



by a; = 2E^/m^. This gives the FSR energy spectrum in 
the rest frame of each for the decay — > //. Next, we 
boost to the rest frame of x, which essentially coincides 
with the Earth's rest frame. The FSR photon spectrum 
in the x I'cst frame is given by 
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which may explicitly be written as 
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Here y = to^ > 2m/, and the upper limit on 

the integral in eq. ^ corresponds to the maximum 
photon energy in the coUinear limit, which approaches 
unity for mf <^ m^. For the decay channel x ~^ 
4'4> ffff: we have instead y = In eq. H]), 

Emax — (-^ — z^r-)my. for annihilating dark matter, and 
the corresponding formula for decaying dark matter has 



Emax = i ™"* ^™^ )^. Note that the factor of two ap- 
pears in eq. ^ since there are two (/)'s produced for each 
dark matter annihilation or decay. 

An important observation is that smaller values of m^ 
lead to a softer photon energy spectrum. We note that 
the FSR formula in eq. ^ receives ) correc- 

tions, which are small in the coUinear limit but which 
become 0(1) when m^ is of the same order as 2m f. We 
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Dwarf 


Vel. disp. 


■^ann 


■^dcc 


ACT 


Flux limit 






[km s~^] 


login [£ann/(GeV^cm-^)] 


logio[£dcc/(GcVcm"^)] 




[cm-^s-'] 




Sagittarius 


11.4 


19.35 ± 1.66 


18.73 ± 1.44 


HESS 


$(S > 250 GeV) < 3.6 x 10" 


-12 


Draco 


10.0 


18.63 ±0.60 


17.51 ±0.12 


VERITAS 


> 200 GeV) < 2.4 x 10~ 


-12 


Ursa Minor 


9.3 


18.79 ±1.26 


17.55 ±0.36 


VERITAS 


> 200 GeV) < 2.4 x 10" 


-12 


Willman 1 


4.3 


19.55 ±0.98 


17.51 ±0.84 


VERITAS 


$(£ > 200 GeV) < 2.4 x 10" 


-12 



TABLE I: Line-of-sight stellar velocity dispersions and integrals over the mass density and density squared for four Milky Way 
dwarf galaxies. Also given are the gamma-ray flux upper bounds for each dwarf from observations with the two Atmospheric 
Cherenkov Telescopes HESS [7l[ and VERITAS [t^]. jC- values given represent on minus off regions as defined in Section UlTl 
and the error bars represent 90% confidence level upper and lower bounds. 



will ignore these corrections, since the only case consid- 
ered in this paper in which they would be important is 
when the dark matter decays or annihilates first into a 
250 McV (f>, which in turn decays into muons. In this 
case, however, the most important contribution to the 
photon energy spectrum is not from FSR off the muon, 
but rather from the radiative decay of the muon, which 
we discuss below. 

Note that we assume there is no contribution from in- 
ternal bremsstrahlung, which we differentiate from FSR 
and define as photons directly radiated from the media- 
tor(s) of the interaction. This is because we assume the 
mediator (s) of both the annihilation and decay channels 
are electrically neutral. 

For the channels that produce muons, we also include 
the effects of radiative muon decay, in which we have the 
decay channels fj,~ — > e~v^Ve'^ and e'^u^Uej. The 

photon spectrum in the muon rest frame has been calcu- 
lated in (tgI in the limit r = <JC 1. For unpolarized 
muons. it is given by (see also [sij ) 
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where w = por the annihilation or decay channel 

x(x) ~^ /^^M J the energy spectrum from radiative muon 
decay in the x rest frame is to a good approximation 
obtained from eq. ([7]) as 
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where y = in the case of annihilations, and y = in 

^ nix 

the case of decays. The factor of two appears in eq. ([8]) 
since two muons are produced in each annihilation or 
decay. 

For the annihilation or decay channel via an interme- 
diate step involving (f>, i.e. xix) ~^ 4>4> ~^ 



two boosts are required. The first boost takes us from 
the rest frame of the muon to the rest frame of the ch: 



where x 
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and B = . 1 #- is not necessarily 

close to one, and a second boost takes us to the rest 
frame of the x- 

2 / dx r^, (10) 



dN^ 



dy 



where again y = in the case of decays, or y = ^ in 
the case of annihilations. Note that for the second boost 
we assumed <C ni^. 

The contribution to the photon spectrum from the ra- 
diative muon decay is subdominant to the FSR off the 
muon for ^ m^, if x goes directly to muons, and 
for 777.0 ^ 777^, if X goes to muons via 0's. For example, 
for XX ^ /^^M I the number of photons above 200 GeV 
coming from radiative muon decay is about (5%, 10%, 
20%) compared to the number coming from FSR off the 



X dx 



muon for 



(500, 10^ 10**) GeV, respectively. (We 



chose to consider photons above 200 GeV since this is an 
appropriate threshold energy for VERITAS.) However, 
the radiative muon decay dominates over the FSR off the 
muon if ^ [3l|, , and also produces many pho- 
tons that are much harder than the FSR contribution. 
For example, for = 250 MeV, the number of photons 
above 200 GeV from radiative muon decay is about 35 
and 10 times larger for 777^ = 10'^ GeV and 10^ GeV, 
respectively, than the number coming from FSR off the 
muon calculated with eq. For 777^ lighter than about 
725 GeV, for the above case, there arc no FSR photons 
off the muon above 200 GeV. However, there are photons 
above 200 GeV from radiative muon decay for even lower 



values of 



Thus, the above-mentioned 0(1) correc- 



tion to the FSR formula in eq. (|3]), which formally may 
be necessary when 777^ ~ 777^, is negligible in this case, 
since the radiative muon decay dominates the signal. 

For the channels that produce r's, we only include the 
cases where x annihilates or decays directly into t~^t~ . 
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In addition to the FSR component, t decays produce 
many tt^'s which each decay into two photons. We take 
the parametrization of the total photon spectrum from 
fr^ . which, for = 1 TeV, these authors found to be 
from simulations 



dy 



-1.31 



(6.94y - 4.93?;^ - 0.51y^)e-'^ '^^y . (11) 



This parametrization is valid down to about y ~ 0.01 
fn\ . where y is defined as above. Although this is the 
parametrization of the spectrum for = 1 TeV, to a 
good approximation the shape of the spectrum as a func- 
tion of y does not change much over the whole range of 
values we consider for (a few hundred GeV to 10 
TeV). In particular, the dominant contribution to the 
photon spectrum coming from the tt*^ decays remains vir- 
tually unchanged, while the sub-dominant contribution 
from the FSR component changes only logarithmically 
as a function of m^, see eq. Note that we also ig- 
nore any uncertainties in extracting the photon spectrum 
from simulations (see e.g. pof). 



III. GAMMA-RAY OBSERVATIONS FROM 
ATMOSPHERIC CHERENKOV TELESCOPES 

In Table |T1 we list the current upper bounds on the 
gamma-ray flux measured for each dwarf galaxy from 
various ACTs. The flux bound for the Sagittarius dwarf 



is from HESS using 11 hours of observation [71|. The 
bounds for Willman 1, Ursa Minor, and Draco are from 
recent VERITAS observations VERITAS observed 
Willman 1 for 15 hours and Draco and Ursa Minor for 20 
hours each, and found gamma-ray flux upper limits for all 
three dwarfs at the level of 1% of the Crab Nebula [t^ . 
We determine the gamma-ray flux of the Crab Nebula 
above 200 GeV, which is the reported energy threshold 
of VERITAS, by integrating the power-law fit for the 
Crab Nebula of 3.2 x lO~^^{E/TeY)-^-^^cm-'^s-^TeY^^ 
found by [73. We list 1% of this value in Table U as 
the flux bounds for Willman 1, Ursa Minor, and Draco. 
Similar flux upper bounds for Willman 1 were foun d by 
the MAGIC telescope over 15.5 hours of observation [79l |. 
MAGIC also observed Draco over 7.8 hours jsO], and 
STACEE observed Draco for a total of 10.2 hours [H. 
Both obtained flux bounds for Draco about an order of 
magnitude larger than VERITAS. 

To obtain the quantities £ann and £doc for each galaxy, 
we must subtract a background level from the signal re- 
gion, where both the signal and background regions are 
specific to each ACT. In what follows we refer to the 
background region as the off mode, and the signal re- 
gion as the on mode. For HESS, the on region is defined 
to be within a solid angle of 0.14° around the center of 
Sagittarius, while the off region is defined as an annulus 
0.43 - 0.57° from the center of Sagittarius. For VERI- 
TAS observations, the corresponding on region is within a 
solid angle of 0.15° centered around each dwarf, while the 
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FIG. 1: Annihilation cross-section upper bounds as a function 
of dark matter mass obtained from gamma-ray observations 
of the four Milky Way dwarf galaxies: Willman 1, Sagittarius, 
Ursa Minor, and Draco. These bounds were calculated using 
the 90% confidence level lower limit on each £a,nn value given 
m Table H] The bounds given are for six different channels, 
which are from top curve to bottom curve: XX ~^ 4"^ ~' 
fj,^ H'^ fj,~ for = 250 MeV, XX ^ 4'4> ^ e^e~e^e~ 
for = 100 MeV, xX ^ /'"^M", XX -> e+e~, xX ^ 4'4> 
for the case where decays directly to 77 with a branching 
ratio of 10% and to e^e~ with a branching ratio of 90%, and 
XX ~^ T^T~ . Here is a new mediator particle as discussed 
in e.g. MM- 



off region is a ring with inner/outer radius of 0.4 — 0.5°. 
In principle, one must include the contribution from the 
diffuse Galactic halo. However, given that we arc con- 
sidering satellites with angular extent much less than a 
degree, the contribution from the diffuse Galactic halo 
emission is similar for both the on and off regions, so our 
results do not depend on the normalization of the diffuse 
halo signal. 

Given the on and off regions above, we give the £ann 
and Cdcc values for each galaxy in Table [D More specifi- 
cally, these are the respective maximum likelihood C val- 
ues for the on region minus the off region. In this way, 
we determine the distribution for the best fitting C val- 
ues for each galaxy, fully accounting for correlated error 
bars. We note that for the angular regions considered 
here, the background flux estimated from the off regions 
equals ^ 10% and 50% of the flux within the on regions 
for annihilation and decays, respectively. For each entry 
in TablclH the error bars represent 90% confidence levels. 
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FIG. 2: Decay lifetime lower bounds as a function of dark 
matter mass obtained from gamma-ray observations of the 
four Milky Way dwarf galaxies: Willman 1, Sagittarius, Ursa 
Minor, and Draco. These bounds were calculated using the 
90% confidence level lower limit on each £dcc value given in 
TableU] The bounds given are for six different channels, which 
are from bottom curve to top curve: x ~^ 4'4' ~^ M^M"^^/^" 
for = 250 MeV, x ^ H ^ e+e"e+e" for = 100 
MeV, X ~* f^^ J X e^e^, x ~^ 4"^ ~* e^e~ for the case 
where (j) decays directly to 77 with a branching ratio of 10% 
and to e^e~ with a branching ratio of 90%, and x ~^ t'^t~ . 



IV. ANNIHILATION CROSS-SECTION UPPER 
BOUNDS 

Figure [T] shows the dark matter annihilation cross- 
section upper bounds as a function of dark matter mass 
for the four dwarf galaxies for six different annihilation 
channels. We used the 90% confidence level lower limit 
on each £ann value to give very conservative cross-section 
upper bounds. The cross-section upper bounds were cal- 
culated using the FSR formulae given in equations ^ 
and ^ , including radiative muon decay for channels that 
include muons. The weakest bounds are for the XX ^ 
(fx/) 1-1'^ fi'^ channel with ~ 250 MeV. (The 
same channel with = 1 GeV (not shown) has a bound 
that is an order of magnitude stronger at ~ 500 GeV, 
and about 25% stronger at to^ — 10 TeV.) Stronger con- 
straints in general are found for the channels that anni- 
hilate directly into leptons without a (j) mediator, since 
here the photon spectrum is harder. However, a strong 
constraint can be found for the xX ~^ 4"^ channel for the 
case where (/> is a scalar of 0(100 MeV) and mixes with 
the Standard Model Higgs. so that (/> decays to e+e" with 
a branching ratio of about 90%, and directly to 77 with a 



branching ratio of about 10% [T^, 112, ISJ] . The strongest 
constraint comes from annihilations into r's, due to the 
large number of high energy photons produced in t de- 
cays to TT^'s. Note that we do not include any boosts 
from halo substructure here, which may strengthen the 
cross-section upper bounds by about a factor of 20 [s^l . 

Among the four dwarfs we consider, the constraints 
from Willman 1 are the strongest. Compared to Willman 
1, the constraints from Draco, Ursa Minor, and Sagittar- 
ius are weaker by a factor of roughly 3.5, 11, and 15, 
respectively. 

We note that in the literature various values of £ann 
have been used to obtain constraints for Sagittarius, 
e.g. £an„ =i (6.4 X 10l^ 4.9 x 10^^ 1.5 x lO^i) GeV^ cm-^ 
assuming a (large core, Navarro-Frenk- White (NFW) 
[sst . small core) profile, respectively - see e.g. [23ll3ll.[3^. 
Our central value of £ann for Sagittarius is close to the 
value quoted for the NFW profile, while the 90% confi- 
dence level lower limit that we use to calculate the upper 
bounds in Figure [1] is about a factor of 13 less than the 
smallest of the quoted values (large core) . 

If we assume that the PAMELA, ATIC, and PPB- 
BETS signals are from dark matter annihilations, then 
one can calculate the dark matter mass and cross-section 
for each annihilation channel that best fits these data 
sets. To fit the PAMELA data alone, a wide range of 
values give a good fit, so we focus on the ATIC/PPB- 
BETS constraints. For the channel xX ~^ e+e", sev- 
eral authors find a best-fit of ~ 700 GeV and 
(crw) - 6 X 10-24 cm^s-i [ll,|2l,[3l|. This cross-section 
is about one order of magnitude smaller than the con- 
servative upper bounds shown in Figure [1] for Willman 
1. For the channel xX ~^ M"*"/^", the best-fit value of 
(crw) - 5 X 10-23 cm^s-i for 771^-1.5 TeV EHHllI, 
is a factor of ~ 5 smaller than the upper bounds from 
Willman 1. 

If we have a mediating the interaction, then the 
channel xX ~^ 4'4> ~* e^e~e^e~ gives best-fit values 
to the ATIC/PPB-BETS data of ~ 1 TeV and 
{av) - 1 X 10^23 cni^s-i [H, [U, [H . This cross-section 
is a factor of ~ 40 below the Willman 1 bounds. For the 
case where <f> decays to 77 with a 10% branching ratio, 
the best-fit cross-section is only a factor of '--^ 5 below 
that from Willman 1. Note that if Sommerfeld enhance- 
ment does not saturate at the local velocity dispersion 
of the Galactic halo, the lower velocity dispersions found 
in dwarf galaxies could boost the cross-section by a fac- 
tor of I'locai /wdwarf ~ 10 to 20 depending on the dwarf. 
Both Willman 1 and Draco would then disfavor this lat- 
ter channel as an explanation of the PAMELA and ATIC 
signals. For the xX ^ 4>4' ^ M^M /^^M^ channel, the 
best-fit value of {av) ~ 7.5 x lO^^s cm'^s"^ for ~ 2.5 
TeV [H [Ml, is a factor of - 50 below the Willman 1 
bound. Again, if Sommerfeld enhancement is relevant at 
the velocity dispersions of dwarf galaxies, then the best- 
fit cross-section could be boosted by a factor of ^ 10 
to 20 in the dwarf, making it close to the current upper 
bounds. 
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Note that if the annihilation cross-section varies with 
velocity through a Sommerfeld enhancement, then a 
more careful treatment of the expected gamma-ray flux 
would allow for the cross-section to vary as a function of 
position within the dwarf galaxy [s^ . This would entail 
including {av) within the integral over the dark matter 
distribution in cq. We find that using this more care- 
ful analysis gives expected boost factors of the same order 
of magnitude, within the solid angles of interest here, as 
we get from a more naive treatment using the ratio of 
the average local and dwarf velocity dispersions. 

Although channels which annihilate into t's do not give 
a very good fit to the PAMELA/ATIC/PPB-BETS data 
p^ . it may still be consistent to fit the data with, for 
example, a 1, 2, or 3 TcV dark matter particle with a 
cross-section of about (3, 10, 20) x lO"^'' cm^s""'^, respec- 
tively [l^, H^. These best-fit cross-sections are ~ 1 to 
2 times larger than the upper bounds from Willman 1 
(and only ~ 2 to 3 times smaller than the bounds from 
Draco), marginally ruling out dark matter annihilations 
into t's as a possible explanation of the data. 

We emphasize that the 90% upper bound on £ann is 
larger than the 90% lower bound by about one to three 
orders of magnitude, depending on the dwarf. It is thus 
perfectly conceivable that further ACT observations will 
see a gamma-ray signal from at least some of these dwarfs 
if dark matter is responsible for the PAMELA/ATIC 
anomalies. 



V. DECAY LIFETIME LOWER BOUNDS 

In Figure O wc give the lower bounds on the dark 
matter decay lifetime as a function of dark matter mass. 
Again we used the 90% confidence level lower limit for 
each iZdec value to give very conservative bounds. The 
channels are the same as those given in Figure [TJ ex- 
cept that one dark matter particle decays into leptons 
or 0's. For decays, the 90% lower bounds are strongest 
from Draco, with Ursa Minor and Sagittarius providing 
the same constraints within a factor of a few. Willman 1 
has about a factor of 5 weaker constraints than Draco. 

If the PAMELA, ATIC and PPB-BETS signals are 
from dark matter decays, then one can calculate the dark 
matter mass and decay lifetime that best fits the data. 
For the channel XX ~^ e+e~, values of ~ 1.5 TeV 
and r ~ 10^^ sec are found to best fit the data sets 
psi . . This decay lifetime is 2 to 3 orders of magni- 
tude above the lower bounds given by the dwarf galaxies, 
depending on the dwarf. For the xX ~^ U~ channel, 



best- fit values of ^ 2 TcV and t ^ 10^^ arc found 
[43! . [46j , again 2 to 3 orders of magnitude above the life- 
time lower bounds given by the dwarfs. Best-fit values 
for yy t^t^ are ~ 5 TeV and t ~ 5 x 10^^ 
{43! . l46j , where the lifetime is about a factor of 20 above 
the current bounds from Draco. We see that the gamma- 
ray constraints from dwarf galaxies do not currently con- 
strain the parameter region of interest for decays. How- 
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FIG. 3: Predicted flux levels (solid black lines) for the Fermi 
satellite from the dwarf galaxy Segue 1, as a function of dark 
matter mass and annihilation cross-section, for four different 
annihilation channels: xX ~^ e^e~ , xX ~^ /^"'"M") XX ~* 4'4' ~^ 
e^e~e^e~ for m0 = 5 MeV, and xX ^ 4"t' ^ M'''A'~A'^M~ for 
7710 = 250 MeV. Here we assume Eth = 100 MeV in eq. HJ, 
and with this Eth, after one year of observing, Fermi can de- 
tect sources above a flux level of 2.4 x 10"^ cm~^s""^ with Scr 
significance (dashed line). The fluxes were calculated assum- 
ing the mean expected value of £ann. For all four channels, 
the dotted lines indicate the approximate values for the an- 
nihilation cross-section and mass suggested by the PAMELA 
data (the ATIC preferred region is, very roughly, a subset of 
this line). For the channels with an intermediate (j), the dot- 
dashed line is the PAMELA suggested region assuming an 
additional maximum Sommerfeld enhancement to the cross- 
section in Segue 1. The expected fluxes for the other dwarf 
galaxies may be estimated by a simple rescaling with appro- 
priate ratios of £ann. 



ever, the 90% upper bound on £dcc can be an order of 
magnitude larger than the 90% lower bound that we used 
to set these constraints (depending on the dwarf), so that 
the gamma-rays from decays to, for example, r's could 
conceivably be seen from dwarfs. 



VI. PREDICTIONS FOR FERMI 

In this section, we discuss the prospects for the Fermi 
Gamma-ray Space Telescope to detect gamma-ray signals 
from dwarf galaxies, assuming that the dark matter anni- 
hilates into charged leptons. Previous work on the Fermi 
detection prospects of a gamma-ray signal from dark 
matter annihilations or decays include e.g. [H, [13, HI], 
and e.g. [s^, [s^, [9l| in the context of dwarf galaxies. 

A disadvantage of the Fermi Large Area Telescope 
(LAT) over the ACTs is that its effective area is only ~ 
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FIG. 4: Predicted flux levels (solid black lines) for the Fermi 
satellite from the dwarf galaxy Segue 1, as a function of dark 
matter mass and annihilation cross-section, for the four anni- 
hilation channels listed in Figure |31 Here we assume Eth ~ 5 
GeV in eq. fl}, and with this Eth, after one year of observ- 
ing, Fermi can detect sources above a flux level of 1.2 x 10"^" 
cm~^s~^ with 3(7 signiflcance (dashed line). The fluxes were 
calculated assuming the mean expected value of .Cann. The 
dotted and dot-dashed lines are the same as in Figure |31 



1 [921 compared with the effective area of the ACTs 
which is ~ 10^ m^. However, an advantage of the Fermi 
LAT over the ACTs, besides covering the whole sky, is 
that the LAT will detect photons with energies down to 
about 100 MeV [o^, while the ACTs only detect photons 
that arc above around a hundred GeV. Since the photons 
expected from dark matter annihilations are never more 
than the mass of the dark matter, Fermi would in prin- 
ciple be able to detect lower mass dark matter particles 
than the ACTs. Fermi would become especially impor- 
tant compared to the ACTs in the case that only the 
PAMELA data is caused by dark matter and the expla- 
nation of the ATIC/PPB-BETS data is due to something 
else. This is because a dark matter particle with a mass 
of 0(100 GeV) can explain the PAMELA data, whereas 
a mass above around 700 GeV is required to explain the 
combined PAMELA/ ATIC/PPB-BETS data. 

Here, we will present the detection prospects for one 
recently discovered particularly promising dwarf galaxy, 
Segue 1 [i^l. This satellite is at a distance of 23 kpc 
from the Sun at a galactic latitude of ~ 50.4° and has 
an observed half-hght radius of ~ 30 pc [ol]. It has 
a measured 1-dimensional stellar velocity dispersion of 
~ 4.3 km/s [62|. Segue 1 has only 24 stars with measured 
line-of-sight velocities , so that its £ value will be less 
constrained than any of the previous four objects wc have 
considered. Using the same formalism as in Section M 



ann,Scg 



^02o.i7±i.44 GeV^cm-s. 



(12) 



Here we have assumed an on region of 0.2° centered on 
Segue 1, and have not performed a background subtrac- 
tion as above when considering ACT limits. This is ap- 
propriate when discussing predictions for Fermi as there 
is no background subtraction in the same fashion as with 
the ACTs. 

In Figure[31 we show with solid black lines the expected 
gamma-ray flux from Segue 1 as a function of the anni- 
hilation cross-section and dark matter mass for the dark 
matter annihilation channels XX e+e~, XX ~^ 
XX 4>4> ^ e~^e~e~^e~ (with = 5 MeV), and 
XX — > 00 — > fi^ fi^ fi^ fi'' (with = 250 MeV). We took 
Eth = 100 MeV in eq. and used the central value of 
/:a„n,scg = lO^^-^GeV^cm-s in eq. In FigureH we 

show the same channels but with Eth = 5 GeV. For both 
figures wc assumed that Fermi will be able to detect and 
reconstruct photons up to 1 TcV in energy. 

To estimate the ability of Fermi to detect gamma-rays 
from a dwarf galaxy requires a careful analysis of the re- 
sponse function of the LAT to the particular gamma-ray 
signal from the dwarf. Such a careful analysis is beyond 
the scope of this work, and we content ourselves with 
a more simple estimate of Fermi's sensitivity, which is 
possible with a few reasonable assumptions. First, we 
assume that the dwarf is a point source. Next, we as- 
sume that the publicly available point-source sensitivity 
as a function of galactic latitude and threshold energy is 
also applicable in this case, even though it was calculated 
assuming a 1/E^ spectrum at the source [92| (the FSR 
signal in eq. ^ has a dependence more like 1/E). With 
these assumptions, the flux sensitivity for Fermi to detect 
a point source at high galactic latitudes with a 3(t sig- 
nificance after 1 year of observation is about 2.4 x 10~^ 
cm-^s-i for Eth = 100 MeV, and about 1.2 x 10"^° 
cm~^s^^ for Eth = 5 GeV. The flux sensitivity is bet- 
ter at higher threshold energies, at least up to about 5 
GeV, above which the sensitivity remains roughly con- 
stant. The background to calculate these sensitivities 
was assumed to be uniform with an integrated gamma- 
ray flux above 100 MeV given by 1.5 x 10~^ cm~^s~"'^sr~^ 
and with a spectral index of -2.1 We indicate the 

respective sensitivities in Figures |3] and S] with dashed 
lines. 

From Figure m we see that generally the 5 GeV energy 
threshold is able to probe more of the cross-section versus 
mass parameter space than an energy threshold of 100 
MeV. The only exception to this is for very low dark 
matter mass (lO's of GcV) for the xX ~^ 00 channel, 
with — > ^.^ fJ-~ and ^ m^. In this case the higher 
energy threshold removes more of the signal than the 
background. 

An exciting prospect is that Fermi could detect a 
gamma-ray signal from Segue 1 if dark matter annihila- 
tions arc responsible for the PAMELA signal and/or the 
ATIC/PPB-BETS signal. In Figures [3] and H the dotted 
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lines indicate the approximate values for the annihilation 
cross-section and mass suggested by the PAMELA data 
for the different channels (the ATIC preferred region is, 
very roughly, a subset of this line). These regions are 
based on the results found in [l5j . and include an ex- 
trapolation to larger masses. (By including only a line 
instead of an extended region, we ignore the uncertainties 
in fitting to the PAMELA data - see e.g. [11].) For the 
channel XX ~^ e+e~, an 0(100 GeV) dark matter par- 
ticle has an annihilation cross-section suggested by the 
PAMELA data that is large enough to give a gamma-ray 
signal detectable with Fermi, assuming a 5 GeV thresh- 
old energy. The cross-sections required to explain the 
PAMELA data with the channel xX ^ A'^M^ f'l'c about 
a factor of two below the one-year Scr detection sensi- 
tivity for masses of a few hundred GeV, but lie right 
at the detection sensitivity for masses of several TeV. 
With a few years of data, the low mass regions can also 
be detected. If dark matter annihilations directly into 
t's (not shown) are responsible for the PAMELA and/or 
ATIC/PPB-BETS excesses, then Fermi can detect the 
resulting gamma-ray signal with about one year of data 
with higher than 5(t significance. 

For the annihilation channels via 0's, the gamma-ray 
signal tends to be smaller so that longer observations with 
Fermi may be required. However, these channels may re- 
ceive an additional Sommerfeld enhancement to the anni- 
hilation cross-section in Segue 1 (up to a factor of ~ 30), 
compared to the local annihilation cross-section required 
to fit the PAMELA data. We indicate the PAMELA pre- 
ferred regions including the additional maximum Som- 
merfeld enhancement with dot-dashed lines in Figures [3] 
and m With the additional boost in Segue 1, the pre- 
ferred cross-sections may well be large enough for the 
gamma-ray signal to be detectable by Fermi. A subtlety 
in calculating the maximum additional Sommerfeld boost 
in Segue 1 is that the enhancement saturates for veloc- 
ities smaller than about m^/m^ [l3|- Saturation does 
not have to occur for the case xX ~* 4'4' ~^ e+e^e^e^ as 
can be as small few MeV, allowing for a large 
additional Sommerfeld boost at the dwarf for all of 
interest. However, for the case xX —^4>4>^ /^^M 

must be at least a few hundred McV, which sug- 
gests that the Sommerfeld enhancement would already 
saturate locally in our halo for dark matter masses of 
only a few hundred GeV. However, for larger dark mat- 
ter masses, an additional enhancement up to a factor of 
~ 30, may occur in Segue 1 (see Figures 3 and 4). 

We note that the discussion in the previous two para- 
graphs assumed that the C value for Segue 1 is the cen- 
tral value given in eq. (|12[) . Since this value has a large 
uncertainty, the prospects for detection may be either 
more pessimistic or more optimistic than presented here. 
Moreover, note that it is easy to rescale the expected 
gamma-ray flux from Segue 1 to estimate the flux from 
other dwarf galaxies. Ignoring the subtleties of the back- 
ground subtraction, which are small for annihilations, 
the rescaling factor is given by the ratio of a particular 



dwarf's C^nn value in Table U to the £ann value of Segue 
1 in eq. . Since £ann is generally smaller for the other 
dwarf galaxies than for Segue 1, the prospects are reduced 
for Fermi to detect a gamma-ray signal from them. How- 
ever, given the uncertainties in the dwarf's dark matter 
distribution, it is conceivable for Fermi to detect a signal 
from, for example, Sagittarius and Willman 1. Moreover, 
it is also possible that as of yet undiscovered nearby dwarf 
satellites with very large mass-to-light ratios could first 
be detected with Fermi through their gamma-ray signal 
from dark matter annihilations. 

In addition to determining the total flux above some 
energy threshold, Fermi may determine the differential 
photon flux d^/dE (which is proportional to dN/dE) 
as a function of energy. Determining d^/dE would be 
helpful in determining whether a detection is a result of 
dark matter annihilation or is from another source. For 
example, the photon spectrum from FSR in the channel 
XX ~^ ff plotted as E'^dN^/dE is unmistakable, since 
it has an edge at the mass of the dark matter particle. 
A detection of such an edge would provide compelling 
evidence that the signal comes from FSR caused by dark 
matter annihilating into leptons, in addition to providing 
a measurement of the dark matter mass. 

Finally, we note that Segue 1 also presents an excellent 
target for ACTs, which in general are more sensitive than 
Fermi in detecting a gamma-ray signal from dark matter 
annihilations in dwarf galaxies if the dark matter mass 
is above the energy threshold of the ACTs. An ACT able 
to achieve the same kind of sensitivity as VERITAS did 
with observations of Willman 1, Draco, and Ursa Minor, 
could potentially detect a gamma-ray signal from Segue 1 
if the dark matter annihilates into charged leptons. Segue 
1 may be a particularly good target for the MAGIC in- 
strument, so we give the background subtracted C values 
for Segue 1 relevant for MAGIC. These are 

C^nn = 102" O4±i-40 GeV^cm^s (13) 

for annihilations, and 

Cdcc = 10i^-80±o-68 GeVcm-2 (14) 

for decays. These numbers assume an on-region of 0.1°, 
and a background region of the same area centered 0.8° 
away from the center of the dwarf, which are the on and 
off regions applicable to the MAGIC telescope [t^ . 

VII. CONCLUSIONS 

We provided conservative bounds on the dark matter 
annihilation cross-section and decay lifetime from final 
state radiation produced by dark matter annihilation or 
decay into charged leptons, and discussed the implica- 
tions of our constraints for the dark matter explanation 
of the PAMELA and ATIC/PPB-BETS excesses. We 
used the experimental gamma-ray flux upper limits for 
four Milky Way dwarf satellites: HESS observations of 
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Sagittarius and VERITAS observations of Draco, Ursa 
Minor, and Willman 1. We find that for the channels 
where dark matter annihilates either directly into elec- 
trons or muons, or through an intermediate particle (res- 
onance), which we call 0, the cross-section bounds from 
the dwarf galaxies do not currently rule out this form 
of dark matter as an explanation of the PAMELA and 
ATIC/PPB-BETS excesses. However, if the mechanism 
of Sommerfeld enhancement is relevant for the factor of 
10 to 20 lower velocity dispersions in dwarf galaxies 
as compared to the local Galactic halo, then the cross- 
sections of dark matter annihilating through 0's required 
to explain the excesses are very close to the cross-section 
upper bounds from Willman 1. Dark matter annihilation 
directly into r's is also marginally ruled out by Willman 1 
as an explanation of the excesses, and the required cross- 
section is only a factor of a few below the upper bound 
from Draco. Dark matter decays as an explanation of 
the observed excesses are currently not constrained by 
gamma-ray observations of these dwarf galaxies. 

We note that since the uncertainty in the gamma ray 
flux coming from the uncertainty in the dwarf dark mat- 
ter distribution is about 1-3 orders of magnitude, de- 
pending on the dwarf, we have provided rather conserva- 
tive constraints by using the lower value of the expected 
fluxes. However, it is possible that further observations 
by Atmospheric Cherenkov Telescopes of various dwarf 
galaxies could see a gamma-ray signal from dark matter 
annihilations [oj (see also [9ll|). 

We also make predictions of the gamma-ray flux ex- 
pected for the Fermi satellite from the dwarf galaxy Segue 
1, which may provide one of the strongest signals from 
dark matter annihilations in dwarfs, and is thus also 
an excellent target for the Atmospheric Cherenkov Tele- 
scopes. We find that in a sizeable fraction of the parame- 
ter space in which dark matter annihilation into charged 
leptons explains the PAMELA excess, Fermi has good 
prospects for detecting a gamma-ray signal from Segue 1 
after one year of observation. 

Note added: While this paper was considered for 
publication in PRD, the measurement of the cosmic ray 
e+ -I- e~ spectrum from 20 GeV to 1 TcV with the Fermi 



Large Area Telescope appeared [95|. Their results call 
into question the sharp peak observed by ATIC/PPB- 
BETS, but they nevertheless still observe an anomalous 
excess that might be due to dark matter. One can again 
calculate the dark matter mass and cross-section for each 
annihilation channel that best fits both the PAMELA 
and Fermi data sets (we use the results in in the 
following). The channel XX ~* e+e" is not a good fit. 
For the channel xX ~^ /^^M i the best-fit value of (av) ^ 
3 X 10-23 cm^s-i for - 1.4 TeV is a factor of - 8 
smaller than the upper bounds from Willman 1. For 
the channel xX ^ r+r^, the best-fit value of (av) 
2 X 10-22 cni3s-i for ~ 3 TeV is a factor of ~ 2 
larger than the upper bounds from Willman 1 , marginally 
ruling out dark matter annihilations into r's as a possible 
explanation of the data. For the channels xX ^ 4'4' ~* 
e^g-e^e" and xX ~* 00 ^ M^M /^^/^ j the best-fit 
values are similar to those for ATIC/PPB-BETS and the 
constraints are thus similar to the values given in Section 
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